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An approximate solution scheme, similar to the Gutzwiller approximation, is presented for the 
Baeriswyl and the Baeriswyl- Gutzwiller variational wavef unctions. The phase diagram of the one- 
dimensional Hubbard model as a function of interaction strength and particle density is determined. 
For the Baeriswyl wavefunction a metal-insulator transition is found at half-filling, where the metallic 
phase (U < U c ) corresponds to the Hartree-Fock solution, the insulating phase is one with finite 
double occupations arising from bound excitons. This transition can be viewed as the "inverse" of 
the Brinkman-Rice transition. Close to but away from half filling, the U > U c phase displays a finite 
Oh Fermi step, as well as double occupations originating from bound excitons. As the filling is changed 

away from half-filling bound excitons are supressed. For the Baeriswyl-Gutzwiller wavefunction 

w : 

at half-filling a metal-insulator transition between the correlated metallic and excitonic insulating 
state is found. Away from half-filling bound excitons are suppressed quicker than for the Baeriswyl 
wavefunction. 
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I. INTRODUCTION 



Variational studies have contributed greatly to our understanding of strongly correlated systems, described by 
the Hubbard model [1-4] and its extensions. While the last decades saw the development of the dynamical mean- 
field theory [5] and the density matrix renormalization group [6] variational studies still play an important role in 
understanding metal-insulator transitions. In part this is due to their relative simplicity and applicability to large 
systems irrespective of number of dimensions. Two frequently used variational wavefunctions, are the Gutzwiller 
wavefunction [1,2] (GWF) and the Baeriswyl wavefunction [7,8] (BWF). The former is based on suppressing charge 
, O , ■ fluctuations in the noninteracting solution, the latter on projecting out fluctuations in the hoppings from a completely 
projected GWF. The combined use of both projectors (Baeriswyl-Gutzwiller wave function (BGWF) or Gutzwiller- 
Baeriswyl wave function (GBWF)) [9,10] has recently raised the possibility of superconductivity in the two-dimensional 
Hubbard model [11,12]. The idea of using projections based on the kinetic energy or more general operators also 
■ appears in continuous models [13-15]. 

\ The GWF can be solved exactly only in one [16] and in infinite [17] dimensions. In one dimension the exact solution 
of GWF is metallic, in contradiction with the exact result for the Hubbard model [18]. The GWF, however was shown 
to be metallic for all finite dimensions [19,20]. In higher finite dimensions high order diagrammatic treatments [21] 
and quantum Monte Carlo [22] are possible. Only extended versions of the GWF can account for insulating behavior: 
when correlations between doubly occupied sites and empty sites are incorporated (bound excitons) [8,23,24], or when 
t-H . the non-interacting wavefunction from which charge fluctuations are projected out is itself insulating [25,26]. 

For the BWF an exact analytical solution is in general not available. It can be shown [20] that the Drude weight is 
identically zero, hence the BWF is an insulating wavefunction. If the Neel state is assumed to be the wave function 
for infinite interaction then a solution is feasible. [20] For the general case analytical approximations exist [7,10,8], 



and quantum Monte Carlo is also applicable [11,12]. In the limit of large interaction it is known that the BWF 
corresponds to bound excitons and is therefore insulating [8]. 

The GWF is often treated via a combinatorical approximation also due to Gutzwiller [2,27-30] (GA). The GA 
predicts a metal-insulator transition (Brinkman-Rice transition) [2,27-30] and is exact in infinite dimensions [17]. 
The relation of the GA to the exact GWF has also been studied [31,32]. In recent work the author and co-workers 
suggested that the GA consists of using a simplified form for the spin correlations in the non-interacting reference 
wave function [32]. Similar approaches, in which the exchange interaction is implemented in an effective way, have 
also been used in continuous systems to obtain approximations for correlation. [33-36] Extensions of the GA include 
the time-dependent case [37] , implementation for the multi-band models [38] , ensembles with varying particle number 
(BCS wavefunction) [39], and the calculation of matrix elements between ground and excited states [40]. The GA 
has also recently been applied to fermions in optical lattices [41-43], combination with DFT (LDA) [44-46], and 
RPA [47-50]. An improved version of the GA was also recently proposed [51]. In the context of high tempereature 
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superconductivity variants of the approximate solution have been applied [52-54] in the resonating valence bond 
(RVB) method [55-57], which is based on a completely projected Gutzwiller wavefunction. 

For the BWF or its extensions (BGWF or GBWF) there has not been an approximate solution of a similar 
type. The aim of this work is to develop such a scheme for the BWF and the BGWF. The assumptions in the GA 
for the spatial distribution are applied here in momentum space, an approximation for the k-space analog of the 
exchange hole, defined as ( v E , G(oo)|ftk<T?T-k'CT' l^ai 00 )), (where ^g(oo) denotes the completely projected Gutzwiller 
wavefunction) is invoked. The phase diagram of the one-dimensional Hubbard model is calculated. The phase 
transions are characterized by a decrease of the Fermi step. At half-filling it is a metal-insulator transition (the Fermi 
step disappears), from an uncorrelated metal to an insulator with finite double occupations. Away from half-filling 
the transition is between an uncorrelated metal (fully localized in momentum space) to a correlated metallic state. 
The correlated metallic state has a finite Fermi step smaller than the Hartree-Fock solution. The double occupation 
tends to the value of the completely projected GWF, but double occupations due to exciton binding are present at 
finite interaction. 

Excitons in metals are rare due to screening by free charge carriers. Recently evidence [58,59] was found for the 
presence of bound excitons in single-walled carbon nanotubes. A carbon nanotube can be seen as a system of low 
dimensionality, hence screening can be expected to be significantly reduced, and the effects of correlations are more 
pronounced . In these systems the experimental absorption lineshape can not be reproduced by a tight-binding 
model alone, many electron effects are included via GW-type approaches [58]. The variational ansatz presented here 
incorporates bound excitons via an approximate variational theory. 

This paper is organized as follows. In the following section the method is presented. In particular the GA in its 
original form is used as a starting point to construct a similar approximation for the BWF and the BGWF. In section 
III the results are presented. Subsequently conclusions are drawn. 



II. METHOD 



A. The Hubbard Hamiltonian and variational wavefunctions 



In this study the Hubbard Hamiltonian [1-4] for spin-unpolarized systems at various fillings will be investigated. 
The Hamiltonian in one dimension can be written 

N N 

H = -t^2 c\ a c ja + u^2 mt n n- (!) 

(i,j)a »=1 

We will assume a system with L lattice sites and with and iVj, particles with spins up and down respectively. The 
idea of the BWF [7,8] is to act with a kinetic energy projection operator on the completely projected GWF. The 
GWF is obtained by projecting out double occupations from a Fermi sea, 



|*g(7)> = cx P (^-7j2 n ^ n ^ \ FS ^ 



(2) 



where \FS) indicates a Fermi sea of non-interacting fermions. 
The BWF can be defined using Eq. (2) as 



JV 



|*u(a)> = cxp ( -a ( ]T c\ a c ]a - ^ CT rW J ] |*c(7 -> «>)). (3) 



Diagonalizing the hopping operator one can also write 



|tf B (a)) = cxp |-aj](e(k) - Av)n k ^ |*g(7 -> °°)>, 



(4) 



with e(k) = — 2fcos(k) and fi a being the chemical potential. The completely projected GWF at half-filling contains no 
double occupations. For finite a, however, double occupations arise as a result of the binding of neighboring up-spin 
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and down-spin particles and their second order hopping processes, as shown by Baeriswyl [8]. In particular Baeriswyl 
has shown [8] that the polarization fluctuations at half-filling have the form 



<* 2 > = 2 E Tp (\ (*g(oo)|S, ■ S,|* G (oo))) , (5) 

(id) V J 

where X denotes the total position operator denotes the spin vector of site i. This expression corresponds to bound 
pairs of double occupations and holes or dipoles with random orientations. Double occupations arise as a result of 
second-order hopping processes. It is also interesting to note that the BWF is closely related to the RVB [55-57]. In 
the RVB a completely projected GWF is acted on by a unitary operator whose exponent consists of a sum of selective 
hopping processes (increase or decrease of double occupations). The approximate solution of this method leads to 
solving a spin-i Heisenberg Hamiltonian, which is also true for the BWF [8]. In the BWF, however, all hoppings are 
included in the projection, hence away from half-filling the charge carriers can be expected to be more mobile. 

The two projections detailed above can also be applied in sequence. Two other variational wavefunctions can be 
obtained by 

l*BG(a,7)> =exp ( E C U'" -J2^ n ™ ) ) I*g(7)>, (6) 

with 7 finite, and 

(7) 



|*Gs(a,7)) = cxp ^-7^ ni t n 4^ 



Eq. (6)(Eq. (7)) is known as the Baeriswyl-Gutzwiller [9,10](Gutzwiller-Baeriswyl) wavefunction. Below, in addition 
to the BWF, an approximation scheme is also developed for the Baeriswyl-Gutzwiller wavefunction (BGWF). 



B. The Gutzwiller approximation 

In the following the essential features of the GA will be given, for details see Refs. [2,27-29]. The GA was developed 
to simplify the sum over determinants that arise when expectation values are evaluated over ^q. The approximations 
are based on the U = 0( 7 = 0) solution. In the position representation the normalization of the GWF can be written 

<* G |*g) =L-^+^)^| 2 )[k;g 7 ]| 2 |D[l;h / ]| 2 exp[-27£>(g / ,h 7 )], (8) 
1 

where the sum is over all configurations of coordinates, g/ and h/ denote the coordinates corresponding to configuration 
I, D(g/,hj) denotes the number of double occupations for the particular configuration, and ^[kjgj] denotes the 
determinant formed of the plane- waves with wavevectors k at positions g/. Due to the determinants only those 
configurations contribute which include up to one particle of a particular spin at each site. One can define the 
unnormalizcd probability distribution in position space, 

/W(g,h) = |2)[k;g]| 2 |S[l;h]| 2 cxp[-2 7 ^(g,h)]. (9) 

Using Eq. (9) one can write relevant expectation values. For example the average double occupation can be written 

To arrive at the GA one replaces the square of the determinants with their averages in the Fermi sea. Considering 
only the up-spin channel one can write the normalization of the Fermi sea as 

t (FS\FS) t = L- N ^\V[k ]gl ]\ 2 = h (11) 
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since the wavefunctions that enter are normalized planewaves themselves. As the sum in Eq. (11) is over all configu- 
rations of up-spin particles on the lattice, such that at most one particle occupies a particular site we can approximate 
each term by its average as 

|S)[k;g / ]| 2 «(|S)[k;g / ]| 2 ) = ^ ) (12) 

where Cjy denotes the number of ways particles can be placed on L lattice sites. The down-spin particles can be 
handled similarly. Substituting Eq. (12) one can write the unnormalized probability distribution in real space as 

P GA (gh)=exp[-2 7 D(g,h)]. (13) 

The average number of double occupations in terms of Pga can be written 

E 7 PGA(g/,h 7 )£>(g / ,h J ) 



, ^ (14) 



but here, contrary to Eq. (10) a constraint has to be introduced over the summations. Only those configurations are 
summed over, which have zero or one particle of a particular spin at each lattice site. 

The approximation for the kinetic energy involves averaging a product of unequal determinants [2,27-29,22,32] 
since the hopping is not diagonal in the coordinate representation. Similar to Eq. (12) this is done by evaluating the 
hopping energy for the Fermi sea, 

T= t (F5|ci t c jt |^) t = L- Jv t^'S)*[k;g' 7 ]2)[k;g 7 ] ) (15) 

where g/ and g'j denote two configurations which differ only in their occupations of site i and j. For g/(g'/) site i is 
unoccupied(occupied) and site j is occupied(unoccupied). The prime on the sum indicates the restriction that gj are 
configurations such that site i is occupied and site j unoccupied. Of configurations with a given pair of sites which 
have one occupied and one unoccupied site there are C^~\. The product of determinants can be approximated by 
the Fermi sea average, since the hopping energy can be evaluated exactly, i.e. 

T=i^*exp[ik-(R i -R j ))]. (16) 
L k 

where Ri and Rj denote the pair of lattice sites involved in the hopping, and the asterisk indicates that the sum be 
performed over occupied states only. The approximation 

2)*[k; §',]!%; g/ ] « <£* [k; g'^k; g/ ]) = T-j-^ (17) 



C 



L-2 
Nf-1 



can be introduced. Using this approximation the average hopping of an up-spin particle from site j to site i can be 
written 



<*l4i- c it!*> C K E'i Pga(Ei, h J )exp[- 7 AD(g / J ,g I ; h 7 

~iL-2 



In Eq. (18) AD indicates the change in number of double occupations due to the hopping. 

Notice that one could arrive at the approximate expressions in Eqs. (13) and (18) via different reasoning [32]. One 
can define a probability distribution over configurations with up to one particle of a given spin at each site and weigh 
each configuration with the weighing factor P ( 3 J 4(g/,h/) = exp(— 27_D(g 7 , hi)). One can define an estimator for the 
hopping energy of the form Texp(— -fAD(g'j, gj; hi)), considering that the hopping operator connects states with 
different number of double occupations. The scaling factor T can be obtained by requiring that the hopping energy at 
U = (7 = 0) corresponds to the kinetic energy of the non-interacting system (Eq. (16)). Below the approximation 
scheme for the BWF will follow these steps. 

The Gutzwiller approximation gives rise to a very simple form for the momentum distribution 

(nk CT ) 7 =n a (l- 9(7)) + g(7)9(/i (T - e k ), (19) 
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where [i a and 0(x) denotes the chemical potential and the Haviside step function respectively, and 



Ej Pga(si, h J )exp[~ 7 A£>(g , J , g7 ; h 7 )] 

9(7H ^ et^^m • (20) 

From Eq. (19) one sees that the momentum distribution at any filling is a constant function with a discontinuity at 
the value of the chemical potential. For half-filling in the limit 7 — > 00 the distribution becomes | for any k, and q(j) 
can be simplified [29] to 



C. Application to the Baeriswyl wavefunction 

The BWF consists of a projection of the fully projected GWF. While the normalization for the GWF can be easily 
written, since the U = solution is known, this is more difficult for the BWF where the U = 00 solution is needed. 
In general one can write the normalization as 

<* B |* B > = 5>(k/,l/)exp -2a^{e(k 7 ) + e(l 7 )} , (22) 
1 Li 

where 

X (kj, 1/) = (* G (7 -»■ oo)|kj, l 7 )(k 7 , 1j|*g( 7 -»■ 00)). (23) 

x(k/, 1/) denotes a probability distribution for a particular set of vectors k 7 and 1 7 which guarantees that at a = the 
7 = 00 momentum distribution (Eq. (19)) is recovered. To account for this distribution one introduces a piecewise 
constant potential discontinuous at \i a in momentum space, so that the distribution in Eq. (19) is reproduced. As 
in the GA the summation in Eq. (22) is such that no two particles of the same spin can occupy the same site in 
momentum space (to account for the Pauli principle), but the distribution is otherwise uncorrelated. The kinetic 
energy is obtained the usual way 



with Pga-b defined as 



P G a-b{^i, h) = X(k/, lj)exp[-2a ]T (e(kj) + e{h))\. (25) 

k 7 ,lj 



In order to arrive at an approximation scheme for the interaction we first write the number of double occupations 
in k-space as 

J2 ni ^ = J (S ft kr fi k';- a kt4<|Ck+ q i-5k'-qi • (26) 

% \k,k' kk'q^O / 

The first term is simply N^N^/L. The second term is a correlated hopping of an up-spin and down-spin particle in 
momentum space. We write 



/ ^ t t \ -E/ J PGA-B(k / ,l / )E kl;ki , q# oexp{-a(e(k 7 )+e(k' 7 )-e(k / +q)-e(k' / -q))} 

\ k £o CktCk ' iCk+qtc " q v = E,^-B(k„i,) 

(27) 

where the double prime indicates that for a particular set ki,kj,q ^ the only configurations which enter are ones 
with k' 7 — q and k'j + q unoccupied and k 7 and k'j occupied. The energy differences in the exponent account for 
the correlated hopping in momentum space. In the original GA applied to the GWF, it is the hopping energy which 
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FIG. 1. Kinetic and potential energies as a function of the variational parameter at half-filling. 




FIG. 2. Kinetic and potential energies as a function of the variational parameter at quarter-filling. 



G 



behaves in a similiar way: there the hopping causes a change in the number of double occupations [2,27-29]. U is 
fixed by requiring that the known number of double occupations is reproduced at U = oo (a = 0) , in other words the 
q = term is cancelled by the q ^ one. Note that in the GA the kinetic energy is multiplied by a scaling factor 
which is fixed by requiring that the noninteracting kinetic energy is reproduced [32]. 

In the limit a — > oo the Hartrcc-Fock non-interacting ground state is obtained. In Eq. (22) the distribution in 
this limit includes only the k vectors corresponding to the lowest eigenvalues. Since the distribution corresponds to 
the finite temperature one with inverse temperature equal to 2a, it also follows that at half filling the discontinuity 
characterizing the Fermi surface of metals is only present when a — >■ oo. For the double occupations, only the q = 
term survives when a — > oo. In the limiting cases U — and U — >• oo the energies are correct within the present 
scheme, as is the case for the standard GA applied to the GWF. Generalization to systems away from half-filling is 
also straightforward, since the distribution for the a — > limit (x(k, 1)) can be chosen accordingly. 



D. Application to the Baeriswyl-Gutzwiller wavefunction 



It is also possible to generalize the above scheme to the combined projection based BGWF [9,10] defined in Eq. (6 
). This generalization is enabled by the fact that the completely projected GWF on which the BWF is based (see Eqs. 
(2) and (3)) enters into the approximation scheme detailed in the previous subsection via the momentum density, 
which for the GA has a known form (Eq. (19)). Hence generalizing the approximation scheme to the BGWF proceeds 
exactly as described above, only that the momentum density is the one corresponding to the value of 7 (Eq. (19)). 



E. Implementation 



In this work our implementation of the above approximations is similar to that described in Chapter 9 of Ref. [29] . 
Assuming that the distribution is fully uncorrelated in k-space one can write the normalization as 

<*b|*b) « lit 1 + (n kCT ) 7 {exp(-2ae(k)) - 1}]. (28) 

kcr 

(n-ka-)'y is the momentum distribution of the Gutzwiller wavefunction (Eq. (19)). The kinetic energy is then obtained 
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k<r kcr 

The resulting interaction energy 



u^imtmjnu-^ (30) 

i 

- ^ exp{-a(e(k) + e(k') - e(k + q) - e(k' - q)}(l - (n k+qt ) 7 )(l - (n k >_ q4 ,) 7 )(n kt ) 7 (n k > i )^ 

where 



nk - k '<' : ' 



C(k, k', q) = (1 + (n kcr ) 7 {cxp(-2ae(k)) - 1})(1 + (n k , CT ) 7 {exp(-2ae(k')) - 1}) (31) 
(1 + (n k+qcr ) 7 {exp(-2ae(k + q)) - 1})(1 + (n k , +q(J ) 7 {cxp(-2ae(k' + q)) - 1}). 

Note that the occupation factors in Eq. (30) are such that states k and k' are occupied, and k + q and k' q are 
unoccupied, which coincide with the correlated hoppings in k-space corresponding to the double occupation operator. 



III. RESULTS 



In Figs. 1 and 2 the kinetic and interaction energies are shown for a twelve site system comparing the results of 
an exact calculation to the outcome of the approach presented here for half and quarter fillings. The kinetic energy 
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shows strong disagreement for intermediate values of the variational parameter, presumably due to the fact that this 
approach does not take into account momentum space correlations. The double occupations are in good agreement 
between the two calculations in both cases. A similar degree of agreement is found at quarter filling. Further testing 
of the method can be seen in Fig. 3 where the energy at half-filling is compared to the exact result [18], the exact 
Gutzwiller result [16,17] and the Gutzwiller approximation applied to the GWF at half-filling. The discontinuity in 
the GA-B results indicates a first-order metal-insulator transition at U c rts 4.04, where the U < U c phase corresponds 
to a — > oo, the Hartree-Fock solution. 




FIG. 3. Comparison of the energy of the GA-B scheme to well known results: GA-G (Gutzwiller approximation applied to 
the GWF), exact Gutzwiller, and exact results. Ubg indicates the transition point for the Baeriswyl-Gutzwiller wavefunction. 



In Fig. 4 the phase diagram is presented calculated using 200 sites. As the density decreases from half-filling, the 
critical interaction strength increases until it reaches a maximum. Similar behaviour is found when the density is 
increased from half-filling. In Fig. 5 the fraction of doubly occupied sites are shown as a function of the interac- 
tion strength at different fillings. For the half-filling case (shown in both panels) double occupations starts at one 
quarter (Hartrcc-Fock value), and then decreases abruptly at U c . Subsequently it decays to zero with increasing 
U. This transition "mirrors" the Brinkman-Rice [27] transition. There, while approaching the critical interaction 
from the metallic side, the number of double occupations decreases. The insulator of the Brinkman-Rice transition 
is the simplest possible insulator, with no double occupations. In the GA-B the double occupations increase when 
approaching the critical interaction from the insulating side, and the metallic side corresponds to the simplest metal; 
the non-interacting Hartree-Fock ground state. 
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FIG. 4. Critical interaction strength as a function of particle density. Closed circles indicate systems in which the large 
interaction phase contains bound excitons. 

Fig. 5 also shows how the fraction of double occupied sites vary for different fillings. For n < 1, close to half-filling 
the double occupations show the same pattern as for half-filling, until at n w 0.8 the state corresponding to large 
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interaction strength no longer contains doubly occupied sites. Above half-filling there is a minimum fraction of doubly 
occupied sites for each system, but close to half-filling we observe a slow tending to the large U limiting value (for 
example n — 1.08). The fraction of doubly occupied sites above the limiting value are due to bound excitons, as 
they would not arise were it not for the Baeriswyl projector. Such bound excitons were only found in the range 
n = 0.8 < n < 1.20. 



0,5 
0,4 

0,2 
0,1 




0.2 



4 


U/t 






5 


n=1.00 


6 


1 


1 1 


n=1.08 








n=1.16 








n=1.20 






1 


















■ 


■ 







0,1 - 



n=1.00 

n=0.92 

n=0.84 

n=0.80 



U/t 



FIG. 5. Double occupation as a function of interaction strength for various fillings. 




FIG. 6. Density as a function of fc-vector at different fillings for various values of the interaction strength. 

In Fig. 6 the momentum density nko- is shown for different fillings as a function of U in the phase corresponding 
to large U in each case. The phases found at small U have a Fermi step of 1 at ir/2. The Fermi step closes entirely 
for the system at half-filling, but remains finite for both systems away from half-filling. Thus away from half-filling, 
we find a a first order phase transition from an uncorrelated metallic phase to a metallic phase which contains bound 
excitons. 

For the Baeriswyl-Gutzwiller projection the approximate scheme presented here results in a minimum energy 
corresponding either to the Gutzwiller or the Baeriswyl wavefunction. For half-filling the transition occurs between 
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a correlated metal and a correlated insulator. The interaction strength at which the transition occurs is given by the 
crossing point of the energy curves GA-G and GA-B, and is indicated in Fig. 3 (Ug G « 6.6). Away from half-filling 
the interaction strength at which the transition occurs increases, and for n = 0.96 we find no transition in the range 
< U < 20: the ground state of the system is a partially projected Gutzwiller function( no bound excitons). For 
n = 0.98 a first order phase transition is found. The transition occurs at Ug G « 7.9. For smaller values of the 
interaction the wavefunction is a partially projected Gutzwiller function, the Baeriswyl projection parameter (a) is 
always zero, only the Gutzwiller parameter (7) varies: the system is a correlated metal without bound excitons. For 
larger values of U the parameter a is finite and approaches zero as U — > 00. The Gutzwiller parameter is such that all 
double occupations are projected out for a = 0, and is constant in this range of U. In Fig. (7) the double occupation 
is shown as a function of the interaction strength. The double occupation decreases linearly for the correlated metal 
described by the Gutzwiller approximation, and is discontinuous at the phase transition. For the large interaction the 
double occupation decays to zero. In the regime where a and the double occupation is finite the double occupations 
can be attributed to bound excitons. Hence away from half-filling a metal-metal transition is found between two 
correlated metallic states, distinguished by the absence or presence of bound excitons. 




FIG. 7. Double occupation as a function of interaction strength for the Baeriswyl-Gutzwiller wavefunction for filling n = 0.98. 



IV. CONCLUSION 



In summary, an approximate scheme was presented to solve the Baeriswyl and Baeriswyl-Gutzwiller variational 
wavefunctions. The approach presented here is simple and easy to apply in finite dimensional systems and large 
system sizes are tractable. The scheme is similar in spirit to the well-known Gutzwiller approximation, in which the 
starting point is the Fermi sea, and the Pauli principle is implemented by requiring that no two particles of the same 
spin can be on the same site in real space, but no other spin correlation effects are included. In the approach described 
herein two particles of the same spin cannot occupy the same site in k-space, hence an approximate treatment of the 
k-space analog of the exchange hole is developed. 

At half-filling a metal-insulator transition is found, where the metallic phase (U < U c ) corresponds to the Hartree- 
Fock solution, the insulating phase is one with finite double occupations corresponding to bound excitons. This 
transition can be viewed as the "inverse" of the Brinkman-Ricc transition. Close to but away from half filling, the 
U > U c phase displays a finite Fermi step (metallic), as well as double occupations originating from bound excitons. 
As the filling is increased or decreased from half-filling bound excitons are supressed. 

For the Baeriswyl-Gutzwiller wavefunction it was found that the optimal solution is always either the Baeriswyl or 
the Gutzwiller wavefunction in this approximate scheme. The phase transitions shift to larger values of the interaction 
strength. At half-filling a metal-insulator transition occurs between a correlated metal (with double occupations 
suppressed) and a correlated insulator (double occupations arising from bound excitons). Away from, but still close 
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to, half-filling a transition was found between two metallic phases, the correlated metallic state arising from the 
Gutzwiller approximation for small interaction, and one containing double occupations arising from exciton binding 
for large interaction. 
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